Coupled dynamics of atoms and radiation pressure driven interferometers 
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We consider the motion of the end mirror of a cavity in whose standing wave mode pattern atoms 
are trapped. The atoms and the light field strongly couple to each other because the atoms form 
a distributed Bragg mirror with a reflectivity that can be fairly high. We analyze how the dipole 
potential in which the atoms move is modified due to this backaction of the atoms. We show that 
the position of the atoms can become bistable. These results are of a more general nature and can 
be applied to any situation where atoms are trapped in an optical lattice inside a cavity and where 
the backaction of the atoms on the light field cannot be neglected. We analyze the dynamics of the 
coupled system in the adiabatic limit where the light field adjusts to the position of the atoms and 
the light field instantaneously and where the atoms move much faster than the mirror. We calculate 
the side band spectrum of the light transmitted through the cavity and show that these spectra can 
be used to detect the coupled motion of the atoms and the mirror. 



I. INTRODUCTION 

In recent years the coupling of the light field inside an 
optical resonator to the mechanical motion of the end 
mirrors 000 has received renewed attention due to 
several developments. The improved capabilities to mi- 
cromachine sub-micrometer sized mechanical structures 
with well defined mechanical and optical properties en- 
ables the creation of moveable mirrors with small damp- 
ing whose motion is strongly affected by the radiation 
pressure due to the light field inside an optical cavity. 
The coupled nonlinear dynamics found in those systems 
is very rich 0, Q and is relevant for systems as large 
as gravitation wave antennas 0, Q> @ down to microma- 
chined cavities 0. The radiation pressure force can be 
used to isolate the mechanical motion of the end mir- 
ror from the seismic background 0. Recently, schemes 
for active cooling of mirrors with radiation pressure have 
been proposed Q and implemented 0, 0, E3- It is 
expected that with these methods cooling of the mir- 
ror motion to its quantum mechanical ground state will 
be possible and that the quantized motion of a macro- 
scopic object, including non-classical states of motion 
and squeezing, can be studied this way [lj, [lj, . Also, 
the mirror could become entangled with the light field 
, and questions related to position measurements near 
the quantum limit can be addressed. In the classical 
regime the mirror motion can exhibit instabilities, self- 
sustained oscillations and chaos due to the retardation 
of the electromagnetic field in certain cases. 

From a theoretical point of view, radiation pressure 
driven cavities offer many challenges and opportunities. 
These systems provide excellent testing grounds for the- 
ories of systems with delay [H . Furthermore quantiz- 
ing the light field is a nontrivial question in principle for 
time-dependent boundary conditions |l8j . 

In this paper we extend the study of radiation pres- 
sure driven interferometers to the case where atoms are 
trapped in the optical dipole potential generated by the 
standing-wave light field inside the cavity, see Fig. 2] In 
this system the atoms, light field and mirror interact with 
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Figure f : (Color online) Schematic of atoms in a cavity with 
a moveable mirror. 



each other on an equal footing in the most general case. 
We derive the equations governing this system from first 
principles. Treating the motion of the mirror as well as 
the light field classically, we show that in many cases of 
interest the motion of the atoms can also be described 
classically. We restrict ourselves to those situations in 
this paper. Our key result is that, surprisingly if one 
considers the complexity of the problem, it is possible 
to develop an effective theory that captures many of the 
aspects of the system and is amenable to a fairly simple 
analysis under moderate assumptions. The key observa- 
tion is that the effect of the atoms on the light field can 
be represented by an effective mirror whose reflection and 
transmission coefficient can be calculated from the self- 
consistently determined atomic density distribution. As 
we will discuss the atoms naturally arrange themselves in 
a configuration corresponding to a Bragg mirror of maxi- 
mum reflectivity for a given number of atoms jl9j |. Thus, 
all degrees of freedom of the atoms collapse into a single 
coordinate describing the position of this fictitous mirror. 

The collective back-action of the atoms, as described 
by what we may call the "atom mirror" , strongly affects 
the light field. Thus the atoms qualitatively change the 
potential in which they move themselves. Furthermore, 
because of the modifications of the radiation pressure 
force, they also influence the mirror motion. 

As a consequence of the back-action on the light field 
the potential wells in which the atoms sit become nar- 
rower and for certain mirror positions exhibit two local 
minima per half optical wavelength, i.e. the atomic posi- 
tion becomes bistable. The character of the bifurcation 
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points where the single potential well branches off into 
two valleys can be modified by changing the cavity pa- 
rameters, such as e.g. the reflectivities of the two end 
mirrors. In that way it is possible to choose which of the 
two valleys is the deeper one and accordingly in which 
direction the atoms will move. The two side valleys can 
also be made equally deep, in which case the branching 
point acts like an atomic beam splitter. We also show 
that the atomic motion can exhibit hysteresis, where the 
atoms move along different paths depending on the di- 
rection of motion of the mirror, giving rise to additional 
instabilities of the mirror motion. 

In this paper we consider only the adiabatic limit where 
the mirror moves slowly compared with the atoms, which 
in turn have a much slower response time than the intra- 
cavity light field. Even in that limit the coupled system 
exhibits a remarkably rich phenomenology that we illus- 
trate by analyzing and discussing the modified dipole po- 
tential felt by the atoms as well as the modified motion 
of the end mirror. 

We propose to measure the sideband spectrum of the 
light transmitted through the cavity as a quantitative 
indicator of the coupled atom-mirror motion, as the cou- 
pling of the mirror motion to the atoms gives rise to 
unambiguous spectral signatures that should be within 
easy experimental reach. 

This article is organized as follows: In Section ITT1 we 
describe the model for the coupled atom-mirror-field sys- 
tem and derive the basic dynamical equations. In section 
IIIII we consider the dynamics of the mirror with field and 
atoms following its motion adiabatically, and discuss sev- 
eral aspects of the atomic motion. In section Hvl we cal- 
culate the spectrum of the light transmitted through the 
cavity and show that it can serve as a convenient probe 
of the coupled mirror-atom motion. Finally, section IVI is 
a conclusion and outlook. 



A. Atoms 

Using Kip as the unit of energy and fi.fco = fi.27r/Ao as 
the unit of momentum the free hamiltonian of the atoms 
takes the form 



Mx), (I) 



where ip and ijfi are the atomic field annihilation and 
creation operators. Our theory can easily be adapted to 
describe fermions or bosons although at the level of our 
present theory there are no significant effects due to the 
statistics of the atoms. For the sake of definiteness we 
nonetheless assume that the field operators in Eq. (JTJ 
satisfy bosonic commutation relations. The energy 
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with M the mass of the atoms is the two-photon recoil 
energy and 



U 



ErA 
X tt ' 



with a the s-wave scattering length is the interaction 
strength between the atoms describing s-wave collisions. 

We assume that the atoms can be approximated as 
two-level atoms with transition frequency tu a and we treat 
their interactions with the light field in the dipole and ro- 
tating wave approximations. We consider only one state 
of polarization of the light field and describe this com- 
ponent by a scalar field. Furthermore we assume that 
the intracavity light field of frequency u>l is far detuned 
from the atomic transition, |A| = \ul — uj a \ flu,"/. 
Here is the Rabi frequency of the transition and 7 
is its linewidth. Then we can adiabatically eliminate the 
upper state and, choosing the z axis along the direction 
of propagation of the light field, the interaction between 
light field and atoms is given by the dipole potential 
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MODEL 



g / <Pxtf(x)1>(x)\E(z)\ 2 \u ± (x,y;z)\ i 



(2) 



The system under consideration is shown schematically 
in Fig. ^ Atoms are trapped in the standing-wave light 
field of a Fabry-Perot cavity with one of its end mirrors 
allowed to move under the effect of radiation pressure 
and subject to a harmonic restoring force, £ being the 
displacement of the mirror from its rest position. The 
cavity length offset L only gives rise to a shift of the 
resonances and we will set it to zero in all that follows. 

We measure all lengths in units of the vacuum wave- 
length of the light and all times in units of fip 1 , where 
ttp is the oscillator frequency of the unperturbed mov- 
ing mirror. The electrical fields are described in terms 
of their intensity integrated over the beam profile, i.e. in 
terms of power. 



where 
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(3) 



where p is the electric dipole moment along the polariza- 
tion direction of the laser field and the electrical field has 
been split into a part varying along the cavity axis z and 
a transverse profile u±(x, y; z) assumed to be normalized 
to unity, 



dxdy\u±(x,y;z)\ 2 = I, 



(4) 



for every z. For simplicity we assume in the following 
that the z-dependence of u± along the atomic cloud can 
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be neglected. We further assume that the light field is 
red detuned, A < 0, so that the atoms are attracted to 
the intensity maxima. 

The standing wave that forms in the cavity results in 
an optical lattice potential. We assume that its wells are 
very deep, i.e. g|-Eo| 2 3> Er where Eq is the amplitude of 
the standing wave. In the language of condensed matter 
physics we assume that the atoms are deeply in the Mott 
insulator regime j2jj. It is then a good approximation 
to assume that every site has a specific number of atoms 
in the ground state of the harmonic oscillator potential 
representing the lattice at that site. Due to the back- 
action of the atoms on the light field the distance between 
the lattice sites is smaller than Ao/2, a point to which we 
return shortly. 

The interaction between atoms only gives rise to a con- 
stant energy shift. Its influence on the localized wave 
functions is irrelevant since all that will be required in 
the following is that the atomic density at every site is 
localized within much better than an optical wavelength. 
The order of magnitude of the localization is given by the 
oscillator length at a single lattice site 
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along the cavity axis and 
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in the transverse direction. For typical laser powers and 
detunings the oscillator length is roughly ten to a hun- 
dred times smaller than the optical wavelength. The os- 
cillator frequencies are 
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for oscillations along the cavity axis and 



g\E \ 2 E R 



(7) 



(8) 



for oscillations in the transverse direction. The param- 
eters that we are interested in are such that the oscilla- 
tor length in the longitudinal direction is much smaller 
than an optical wavelength and the oscillator length in 
the transverse direction is much smaller than the beam 
waist. Since the beam is typically much wider than an 
optical wavelength the atoms can be thought of as form- 
ing a stack of pancake shaped discs. 

We consider an optical lattice with N sites occupied 
by n atoms each. Due to the phase shift suffered by the 
light upon propagation through each atomic sheet the 
lattice sites move closer together, the self-consistent pe- 
riod of the resulting lattice being determined later on. 
If the temperature of the atoms is far below the oscil- 
lator energy, we are justified in neglecting thermal exci- 
tations. If the atoms start in the lattice ground state, 
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Figure 2: Field geometry at movable mirror. 



the atomic motion is then reduced to the motion of a 
Gaussian wavepacket in a slowly changing harmonic po- 
tential, a situation that is well described by the classical 
dynamics of the atomic center of mass. 

Because atoms at different sites see the exact same 
potential at all times they will be at the same displace- 
ment from their respective local minima given that they 
started at the bottom of their respective potential wells. 
Therefore it is possible to describe the atomic motion in 
terms of a single coordinate z a , the position of each atom 
modulo the lattice period. 



B. Moving mirror 

The mirror motion is described by the Newtonian 
equation of motion 



mpQpXo 



acF RP , 



(9) 



where c is the speed of light, T" 1 is the mechanical quality 
factor of the mirror, Frp is the radiation pressure force, 
and 

a^ 1 = nipQpXoc 

is a power characteristic of the mirror. The electromag- 
netic field is treated classically and from Fig. [21 we read 
off that 

Fkp = -(\Ei\ 2 + \E 3 \ 2 -\E 4 \ 2 -\E 2 \ 2 ) 
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where we have used that in our situation En = 0. 



(10) 



C. Light field 

We now turn to the back-action of the atoms on the 
light field. The influence of the cavity mirrors will be 
considered in the next section. 

It is well known [2l| that for large detuning, optically 
driving the atoms results in the polarization 
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leading to the inhomogeneous optical field propagation 
equation 



R, T 



V 2 {E Ul _) + k 2 Q (E Ul _) = X \ip\ 2 (E Ul _) 
where in our units the polarizability is given by 
2/j, n P p 2 uj 2 L 
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In principle, the shape of the transverse profile u± is 
the result of the complicated dynamics resulting from the 
repeated propagation of the light field through the cav- 
ity and its interaction with the complicatedly structured 
dielectric constituted by the atoms. A complete descrip- 
tion of these processes is beyond the scope of this paper. 
We instead adopt the point of view that u± has been es- 
tablished in some way and is assumed to be known. We 
therefore describe it an effective manner and eliminate it 
from the theory. 

Specifically, we assume that u± varies slowly with z, 
\du±/dz\ <C |fcoWj_| so that we can neglect its derivatives 
with respect to z. Then, after dividing by Eu±_ to sepa- 
rate the variables, we have from the wave equation 



Q2E ( Z ) , ,2 
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+ k 2 E{z) = (3 2 (z)E(z) 



(14) 



where /3(z) is determined for every z by the eigenvalue 
problem 

X\tp{x, y, z)\ 2 u±(x, y; z)-V 2 ± u_l(x, y; z) = l3 2 (z)uj_(x, y; z). 

In principle this formula solves our problem since it ex- 
presses the propagation equation for E in terms of ip and 
u± independently of how the latter was obtained. In 
order to obtain a more transparent formula we multiply 
Eq. (|15|) by , integrate over x and y and after a partial 
integration we find 



P 2 {z) = XJ dxdy\u±(x,y,z)\ 2 \ij)(x,y,z)\ 2 + 

+ J dxdy(V±u*j_(x,y;z)) ■ (X7 ± uj_(x, y; z)) 
= xWW+k 2 ± (z), (16) 
ending up with the intuitively appealing equation 

d 2 E(z) 



dz 2 



+ (k 2 -ki(z))E(z)= X mz)\ 2 E(z), (17) 



where li/K 2 )! 2 can be interpreted as the density distribu- 
tion averaged over the transverse beam profile. In the 
following we neglect the reduction of the wave vector k± 
due to the finite extend of the beam in the transverse 
direction, an approximation valid if the beam is not too 

strongly focused. 

For the situation at hand |?/>(z)| 2 is readily found, since 
the wave function of the harmonic oscillator factorizes 
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Figure 3: Schematic of the Bragg mirror formed by the atoms. 
Each atomic sheet has a reflection and transmission coefficient 
R s and T s resulting in a total reflection and transmission co- 
efficient R and T. 



into a z-dependent part and a transverse part. As dis- 
cussed above the atomic density distribution is much nar- 
rower than the beam waist w. Hence we find, assuming 
a Gaussian beam profile, 
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where ip z (z) is the z-dependent factor of the oscillator 
wave function of the atoms. 

Equation fT7| allows us to determine the effect of the 
atoms on the light field. The basic idea is to replace the 
entire configuration of atoms by a fictitious "atom mir- 
ror" at some position z a and hence to describe them in 
terms of a pair of effective reflection and transmission 
coefficients. The period of the self-consistent lattice au- 
tomatically fulfills the Bragg condition. Due to the self 
consistency requirement this is even true for a lattice with 
nonuniform filling, see Ref. 0] for more details. (Par- 
enthetically, this means that we can take into account 
an external trapping potential through an appropriate 
choice of effective N and n. The thermal motion of the 
atoms will lower the reflectivity of the atom mirror and 
could be taken into account through a Debye- Waller fac- 
tor.) All that remains to do to describe the effect of the 
atoms on the light field is to determine the transmission 
and reflection coefficient T and R of the atom mirror. A 
schematic of the situation is shown in Fig. [3] 

Consider first a single atomic sheet. Making use of 
the fact that the atoms are localized in the z-direction to 
a much narrower region than an optical wavelength, we 
readily find the transmission and reflection coefficient 
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(19) 



(20) 



by matching boundary conditions at the atom sheet. 
Here 



arctan A 



(21) 



is the phase shift suffered by light upon transmission 
through the sheet of atoms and we have introduced the 



5 



parameter 



A = A (22) 

2TTk Q W 2 K ' 

for notational convenience. Equation (|21|l allows us to 
determine the self-consistent period of the optical lattice. 
Due to the interaction between the light field and the 
atoms the period of the optical lattice is reduced to 

7T - 2(j> 



d = 



(23) 



The entire lattice of atoms can be considered as a peri- 
odically stratified dielectric medium. Its total reflection 
and transmission coefficients are readily found using stan- 
dard methods of the optics of dielectric films |l9|, |22j . In 
essence, one can find the transfer matrix of a unit cell of 
the lattice from the transmission and reflection coefficient 
Eq. H2UI19|1 • The transfer matrix of the entire lattice is 
then the A^th power of the transfer matrix of the unit 
cell. This can be done in closed form using Tchebychev 
polynomials since we neglect scattering losses and there- 
fore the transfer matrix of the unit cell is unimodular. 
From the total transfer matrix the total reflection and 
transmission coefficient are extracted as 251 



R 



T = 



-iANe 2 ^ 
1 - iAN 

p 2i<j>N 



(24) 
(25) 



1 — iAN 

An example of the transmission and reflection coefficient 
of such a lattice of atoms is shown in Fig. 0|for the case of 

Rb. From the figure it is apparent that rather high re- 
flectivities can be obtained in principle. One of the main 
experimental difficulties will be to achieve fairly high fill- 
ing factors over many sites. Also, the high reflectivities 
in Fig. 01 were obtained by using a small beam waist 
and a fairly small detuning of around 100 lincwidths. It 
should be noted that already fairly moderate reflection 
coefficients of the atoms of around 10% give rise to inter- 
esting effects as we discuss later in this paper. S. Slama 
et al. have recently reported Bragg reflection coefficients 
of the order of 30% [2ij. (Note that the assumption of 
a constant and small beam waist is not as bad as might 
be expected, as the Gaussian shaped atomic wavepack- 
ets at every site act like lenses that constantly refocus 
the beam. A more detailed analysis of the spatial mode 
structure of the light field inside the cavity will be pre- 
sented in a future publication.) 

In the following we consider the three representative 
cases of an empty cavity, a cavity containing atoms giving 
rise to an intermediate reflection coefficient of a few tens 
of a percent and atoms with a reflection coefficient close 
to one. 
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Figure 4: (Color online) Reflection (blue dashed line) and 
transmission coefficient (red solid line) of an optical lattice 
consisting of ten 87 Rb atoms at every site as a function of 
the number of lattice sites for A = 10 9 s~ 1 and w — 10Ao as 
obtained from Eqs. I24I25H . 








Figure 5: Representation of the atoms as a mirror at z a and 
illustration of the steady state boundary conditions at z a . 



motion of the atoms which is in turn very slow compared 
to the dynamics of the light field. For this to be true we 
must have Up <C w osc z and w OSCiZ <C k where n is the 
linewidth of the cavity. This is the most relevant case for 
the typical systems currently available. The light field 
then follows the atoms adiabatically, so that both light 
field and atoms adjust instantaneously to any position 
of the moving mirror. We can then eliminate the atoms 
and the light field and express the radiation pressure force 
acting on the mirror in terms of the mirror position alone. 

To obtain the effective radiation pressure force in pres- 
ence of the atoms we observe that, given positions of the 
mirror and of the atoms, the equilibrium light field is 
determined by the boundary conditions at the fictitious 
"atom mirror," see Fig. [SJ The steady-state boundary 
conditions are found to be 
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ADIABATIC DYNAMICS 
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We restrict ourselves to the limit where the motion 
of the moving mirror is extremely slow compared to the 



from which we obtain, e.g. for £4, 
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From these fields we find the adiabatic dipole potential felt by the atoms for a given displacement of the moving 
mirror as 



V{z a -0 = g\ti + ^ 
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a result that holds for every lattice site provided that 
Ri ~ R 2 . Otherwise the depth of the potential wells is 
different at each lattice site, pfij 

The potential V{z a ;£) is periodic with period Ao/2 in 
both z a and £ direction. In the numerator we have, for 
I-R2I close to one, the well known cos 2 k (L + £ — z)- 
type potential. Its minima are moving as the cavity end- 
mirror is moving. 

The effect of the denominator, which gives rise to a 
series of resonances, is somewhat less intuitive. If the re- 
flection coefficient of the atoms is small their back-action 
on the cavity field is weak and we see essentially the res- 
onances of an empty Fabry-Perot cavity: Whenever the 
distance between the two cavity mirrors is (m + ^)Ao/2 
with m an integer the field inside the cavity is very high 
and accordingly the dipole potential very deep [27| ■ Fig. 
shows how a small atomic reflectivity slightly distorts 
this type of potential. The zeros of the numerator are 
visible as dark diagonal lines and the nonzero regions are 
structured by the back-action of the atoms on the light- 
field. 

The situation is qualitatively different if the reflection 
coefficient of the atoms is high, see Fig [7| In this case 
the dipole potential is deep whenever the atoms are in 
a position where they form a resonant cavity with one 
of the end mirrors. Interestingly, this typically gives rise 
to two potential minima per Ao/2, resulting in two dis- 
tinct modes of motion for the atoms: They can follow 
the moving mirror in the local minimum along diagonals 
in Fig. [7| or alternatively they can remain locked with 
the fixed mirror in the local minimum following vertical 
lines in Fig. \7\ These two minima arc actually already 
present for small atomic reflectivity, but only for a very 
narrow range of £ values near £ « 0.25. 

The relative depth of the two local minima is deter- 
mined by the reflectivities of the two end mirrors. The 
potential well closer to the mirror of higher reflectivity is 
deeper. As illustrated in Fig. [5] the system shows behav- 
ior reminiscent of an "avoided crossing" where the two 
resonances corresponding to the two modes of motion 
cross: In positions where the atoms would be on reso- 
nance with both mirrors the potential depth has a fairly 
narrow local maximum. The width of the potential bar- 
rier arising this way becomes smaller as the reflectivity of 
the atomic mirror increases. If the moving mirror sweeps 




Figure 6: (Color online) Dipole potential felt by the atoms for 
Ti = 0.1, T-2 — 0.2 and relatively high transmissivity of the 
atomic mirror, T = 0.90. Here, and in the potential plots to 
follow, lighter regions are deeper in potential and the contour 
lines are drawn at logarithmic intervals for clarity and, as in 
the rest of this paper, all lengths are given in units of the 
wavelength of the injected light field. The green dashed line 
illustrates the path the atoms take when the mirror is moving 
starting at the green dot with the label 1 while the blue solid 
line shows the path the atoms take when the mirror moves 
back down starting at the blue dot labeled 2. 



through such a resonance, the atoms have to decide which 
branch of the potential valley to follow on the other side 
of the resonance [28j. Since we can decide which branch 
of the potential well is deeper by choosing the reflectiv- 
ities of the end mirrors appropriately we can select the 
mode of motion that the atoms will perform. 

The adiabatic motion of the atoms is particularly inter- 
esting if the fixed mirror has the higher reflectivity. We 
discuss it in detail for the case of small atomic reflectiv- 
ity. The trajectories of the local minima corresponding 
to the mirror moving in the positive, resp. negative £ di- 
rection are drawn as a green dashed line and a blue solid 
line, respectively, in Fig. [S] If the moving mirror starts 
in a region far off resonance, e.g. at the green dot labeled 
"1" in Fig. there is only one potential minimum corre- 
sponding to the maximum of the numerator in the dipole 
potential. Moving in the direction of growing £ the po- 
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Figure 7: Dipole potential for higher atomic reflectivity corre- 
sponding to T = 0.05 and Ti = 0.2, T 2 = 0.1. For this choice 
of reflectivities of the cavity end mirrors the potential well 
closer to the moving mirror is always deeper and, in contrast 
to Fig. |S]the atoms move on diagonals following the bottom 
of the potential well for both directions of the mirror motion. 



Figure 9: (Color Online) Radiation pressure force for the two 
directions of motion of the mirror. The atoms take the paths 
shown in Fig. |S| for the mirror moving in the positive and 
negative ^-direction, respectively. The difference in area of 
the two curves is a measure of the work performed on the 
moving mirror. All parameters as in Fig. [(J 



0.05 



0.04 



0.03 



0.02 



0.01 



-0.01 



-0.02 



-0.03 



-0.04 F 



-0.05 




0.24 0.26 0.28 0.3 



Figure 8: Enlargement of the resonance region of Fig. [7| 
showing that the potential well corresponding to motion with 
the moving mirror is deeper than the one corresponding to 
the resonance with the fixed mirror, and explaining why the 
atomic motion is locked to the mirror motion, the figure also 
shows the avoided crossing of the dipole potential. 



tential valley eventually reaches a branching point and 
the atoms follow it to the left towards the fixed mirror, 
since the potential well to the left is deeper if the fixed 
mirror has the higher reflectivity. Near the diagonal zero 
lines of the numerator the potential well becomes shal- 
lower and shallower until finally it ceases to provide a 
local minimum. The atoms "roll" to the second branch 
of the potential valley closer to the moving mirror and 
follow this local minimum until the system hits the next 



resonance, where this story repeats itself. 

The atomic motion is totally different if the moving 
mirror moves in the negative £ direction. Then the 
branch of the potential closer to the moving mirror pro- 
vides a local minimum all the way and the atoms never 
move to the left branch, see the blue line in Fig. 

The motion of the atoms is similar for high atomic 
reflectivity If the mirror is moving in the positive £ di- 
rection the atoms turn left at every resonance into the 
vertical potential valley. Just before reaching the zero 
potential line they roll to the right and fall into the di- 
agonal valley. If the mirror moves back in the negative £ 
direction the atoms also turn left at every resonance and 
move to the left into the diagonal potential valley before 
reaching the zero potential line. 

Interestingly, the radiation pressure force is different 
for the two distinct paths the atoms take as the mirror is 
moving back and forth, see Fig. [5J Hence the light field 
performs work on the mirror in every round trip, and this 
work can be fairly large. It should be noted however that 
our model is too simplistic to study this feature quantita- 
tively. First, the assumption of adiabaticity of the atomic 
motion is likely to break down near the zeros of the nu- 
merator, where the atoms take a sharp turn to the right. 
The potential well becomes extremely shallow in this re- 
gion and when the atoms fall into the other branch of 
the potential, both situations for which the assumption 
of adiabaticity of the atomic motion may not hold. If the 
atoms have some inertia, they could overcome the poten- 
tial wall that prevents them from going over the potential 
hill and continue their motion on the other side. Tunnel- 
ing through the barrier is another possibility. For these 
reasons we concentrate in the following on the case where 
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Figure 10: (Color online) Radiation pressure force for atomic 
motion locked to the motion of the moving mirror for in- 
termediate atomic transmissivity T = 0.9 (red dashed line), 
small atomic transmissivity T = 0.05 (green solid line) and 
an empty cavity (blue dash dotted line) for comparison. To 
ensure that the atoms are following the desired branch of the 
potential we choose T\ = 0.2, T2 = 0.1. The injected laser 
power is 0.01W and a = WW' 1 . 
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Figure 11: (Color online) Radiation pressure potential for 
atoms moving with the moving mirror and an empty cavity. 
All parameters as in Fig. 1101 



of the potential is near £ ss 5 in this case. They are 
smaller for the double peaks for high atomic reflectivity 
and smallest for the empty cavity. 



the atoms move with the moving mirror, i.e. where the 
reflectivity of the moving mirror is higher than that of 
the fixed mirror. For this type of motion the adiabaticity 
assumption holds at all times. 

Once the position of the atom mirror z a (£) has been 
determined we can find the radiation pressure force by in- 
serting z a (£) into the expression for 64. Fig. 1 101 shows the 
three different cases of low atomic reflection coefficient, 
high atomic reflection coefficient as well as an empty cav- 
ity for comparison. A small atomic reflectivity broadens 
the resonance. Also, the peak radiation pressure force 
increases. As the atomic reflectivity increases the reso- 
nance becomes wider and wider and develops a broad dip 
in the middle. For very high atomic reflectivities adjacent 
resonances "collide" and give rise to an avoided crossing 
behavior that results in the double peak structure of Fig. 
I1UI The double resonances in this regime are much wider 
than the empty cavity resonances and the peaks are sub- 
stantially higher. Also the radiation pressure force never 
falls off to values as low as for the empty cavity. The low 
intensity region between the double peaks are the limit- 
ing shape of the aforementioned dips that develop in the 
resonant lineshape. 

From the radiation pressure force we find the total po- 
tential of the moving mirror by a simple integration and 
by adding the harmonic potential due to the restoring 
force. Figure ITTI shows the total potentials, again for the 
three different cases considered for the radiation pressure 
force above. Every time the mirror moves through a res- 
onance the potential makes a steep step down. The steps 
are highest if the atoms have intermediate reflectivity due 
to the broad and high resonances. The global minimum 



IV. SIGNATURES OF THE COUPLED MOTION 

The question of how to detect the coupled motion of 
the mirror and the atoms naturally arises. A possible ob- 
servable is the spectrum of the light transmitted through 
the cavity. The amplitude and phase modulation of the 
optical field due to the mirror motion are detectable as 
sidebands in the spectrum of transmitted light. 

To find the spectrum of the transmitted light we inte- 
grate the Newtonian equations of motion of the mirror 
subject to the modified radiation pressure force using the 
velocity Verlet algorithm. [2p| 



A. Large amplitude oscillations 

First we analyze what happens as the mirror undergoes 
large amplitude oscillations such that in every period it 
sweeps over several resonances. _30j We consider again 
the three cases of atoms with high reflectivity, interme- 
diate reflectivity and an empty cavity. The cavity and 
atomic parameters are the same as before and as for the 
quality factor of the cavity we choose T _1 = 10 4 . 

Figure ^| shows the resulting sideband spectra. For 
clarity each data point has been obtained by averaging 
over 1000 frequency bins. The spectrum of the empty 
cavity falls off nearly exponentially, a consequence of the 
near-Lorentzian lineshape of the resonances and the fact 
that the mirror moves through each resonance at nearly 
constant velocity. 
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Figure 12: (Color online) Sidebands ol the transmitted light 
for atoms moving with the mirror, staying with the fixed mir- 
ror and for an empty cavity. 



The widths of the transmission spectra are compara- 
ble in all three cases, confirming that the widths of the 
resonances are of the same order. However the spectra 
acquire a peak near u> — for nonzero atomic reflec- 
tivities, a consequence of the larger resonance widths in 
these cases; in addition the light field never falls off to 
zero, in contrast to the empty cavity case where it does 
so to a good approximation for high-reflectivity mirrors. 
Especially for high atomic reflectivity the cavity always 
leaks a significant amount of light that gives rise to an 
enhanced DC component in the transmission spectrum. 

The fringes that appear in the case of high atomic re- 
flectivity can be understood in the following way: two 
light pulses exit the cavity whenever the mirror passes 
through a resonance, see Fig. 10, and they give rise to a 
double-slit diffraction pattern in frequency space. Their 
spacing corresponds to the inverse of the time it takes the 
mirror to travel from the first peak to the second peak of 
the double resonance. 

The inset of Fig. ^| shows the details of the low fre- 
quency part of the spectrum for the case of high atomic 
reflectivity. (That low frequency part is similar in all 
cases.) It consists of a peak at the frequency correspond- 
ing to the time it takes the mirror to move from one res- 
onance to the next and of a series of additional peaks at 
each harmonic. These strong components are a result of 
the highly nonlinear potential. The fact that the mirror 
moves at different speeds between different resonances 
gives rise to a broadening of the harmonics into bands 
whose width increases with growing harmonic number to 
finally overlap. Effects due to the damping of the mirror 
motion cannot be resolved because of this "inhomoge- 
neous broadening" . 
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Figure 13: (Color online) Approximate harmonic frequencies 
near the local minima for atoms with high reflectivity, inter- 
mediate reflectivity and an empty cavity. All parameters as 
inFigEl 



B. Small amplitude oscillations 

In our units the average thermal energy of the mir- 
ror at room temperature is T sa 10~ 5 . Comparing with 
the energy scale in Fig. 1111 we see that the mirror can 
undergo many oscillations in each one of the sawtooth- 
shaped local potential minima before it hops to a different 
local minimum due to thermal activation or some other 
perturbation. Thus it is an interesting question to deter- 
mine the oscillation frequency at every local minimum 
for small amplitude oscillations. 

Figure ^| shows these frequencies for the three cases 
already discussed. For high atomic reflectivity, and in 
contrast to the intermediate reflectivity and empty cav- 
ity cases, the oscillation frequency near £ = is much 
lower than the restoring frequency of the mirror. This is 
because the resonances have wide wings in this case and, 
if added to the harmonic oscillator potential of the mir- 
ror, they can compensate the harmonic restoring force 
almost completely over a large range of displacements £ 
and give rise to an extremely flat potential near the local 
minimum. In the other two cases the onset of the radia- 
tion pressure force near the resonances is very abrupt and 
the frequency of the mirror motion can never be signif- 
icantly smaller than the oscillator frequency. At higher 
displacements from the mirror rest positions the stronger 
resonances found with atoms inside the cavity give rise to 
stronger confinement and lead to resonance frequencies 
that eventually become higher than in the empty cavity 
case. Note that for the empty cavity cases we could find 
only a few more local minima beyond £ = 8 and these 
could not be reasonably approximated by a harmonic po- 
tential. 
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V. CONCLUSION AND OUTLOOK 

We have developed a model that describes the coupled 
motion of a cavity end-mirror and atoms trapped in the 
standing-wave intracavity field. The key result is that 
the resulting model is simple enough to allow for a fairly 
straightforward analysis. The self-consistent electric field 
inside the cavity containing the atoms plays a central role 
since from this field we find the modified dipole potential 
in which the atoms move as well as the modified radiation 
pressure force acting on the moving mirror. We find that 
as a result of the strong coupling between light field and 
atoms due to the collective back-action of the atoms the 
dipole potential changes qualitatively, the most notable 
feature being that the atomic position becomes bistable. 
The resonances of the cavity are broadened and develop 
a double peak structure in the limit of high atomic re- 
flectivity. We have studied the dynamics of the mirror 
subject to the modified radiation pressure force in the 
adiabatic limit where the atomic motion is much faster 
than the motion of the mirror and the light field is slaved 
to both atoms and mirror. We have shown that the side- 
bands of the light transmitted through the cavity can 
serve as an indicator of the coupled motion of the mirror 
and the atoms. 

Many open questions remain to be studied. First of 
all a more careful analysis of the spatial mode structure 
inside the cavity in the presence of the atoms is desir- 
able. Such an analysis could start from the Maxwell- 
Bloch equations and should include scattering losses due 
to spontaneous emission from the atomic excited state. 
It would yield more accurate values for the effective re- 
flection coefficient of the atoms and the modified dipole 
potential, especially if R\ and R2 are very different. 

One should also consider different adiabaticity regimes, 
or abandon the assumption of adiabaticity altogether. 
The situation where the atomic motion is roughly as fast 
as the mirror motion and the atom and mirror influence 
each other on an equal footing is particularly interest- 
ing, and so is the case where the retardation of the light 
field is taken into account. Already in an empty cavity 
the retardation of the light field gives rise to many in- 
teresting phenomena such as instabilities, self sustained 
oscillations and chaotic behavior. It will be interesting 



to see how all these effects are influenced by the presence 
of the atoms. 

A related question of practical significance is what hap- 
pens if the atoms are subject to a friction force. Atoms 
can be efficiently cooled to very low temperatures with 
laser cooling and one could hope that by means of the 
coupling to the mirror through the cavity field one should 
be able to also cool the mirror. Loosely speaking, if the 
atoms are subject to a friction force, they lag slightly 
behind the position at which they ought to be accord- 
ing to the position of the mirror. This retardation has 
been seen to lead to the damping of the mirror motion 
in preliminary results, but these calculations need to be 
carried out in much more detail, especially since the mo- 
mentum kicks in laser cooling, and hence the fluctuations 
in the atomic position, are substantial and cannot be ne- 
glected. They can be taken into account e.g. using a 
master equations approach. 

Finally, one needs to address the situation where all 
constituents of our model are quantized. Quantization of 
the atomic motion should allow one to study whether the 
resonances can be used as beam splitters, and the prop- 
erties of the resulting interferometer. Also, a detailed un- 
derstanding of the role of damping of the atomic motion 
requires a proper quantum treatment. The quantization 
of the mirror motion is called for in view of the the exper- 
imental goal to cool it to its quantum mechanical ground 
state. Finally, a quantum treatment of the light field, 
and especially of its fluctuations, is interesting because 
of the powerful methods to measure the fluctuations of 
the electromagnetic field. If they are linked to fluctu- 
ations of the atomic position or the mirror position, as 
must be expected, one could then use the light field as a 
powerful diagnostic tool for the study of these mechanical 
fluctuations. 
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